TASEP WITH DISCONTINUOUS JUMP RATES 



NICOS GEORGIOU, ROHINI KUMAR, AND TIMO SEPPALAINEN 

Abstract. We prove a hydrodynamic limit for the totally asymmetric simple exclusion 
process with spatially inhomogeneous jump rates given by a speed function that may admit 
discontinuities. The limiting density profiles are described with a variational formula. This 
formula enables us to compute explicit density profiles even though we have no information 
about the invariant distributions of the process. In the case of a two-phase flux for which 
a suitable p.d.e. theory has been developed we also observe that the limit profiles are 
entropy solutions of the corresponding scalar conservation law with a discontinuous speed 
function. 



1. Introduction 

This paper studies hydrodynamic limits of totally asymmetric simple exclusion processes 
(TASEPs) with spatially inhomogeneous jump rates given by functions that are allowed to 
have discontinuities. We prove a general hydrodynamic limit and compute some explicit 
solutions, even though information about invariant distributions is not available. The 
results come through a variational formula that takes advantage of the known behavior 
of the homogeneous TASEP. This way we are able to get explicit formulas, even though 
the usual scenario in hydrodynamic limits is that explicit equations and solutions require 
explicit computations of expectations under invariant distributions. Together with explicit 
hydrodynamic profiles we can present explicit limit shapes for the related last-passage 
growth models with spatially inhomogeneous rates. 

The class of particle processes we consider are defined by a positive speed function c(x) 
defined for lower semicontinuous and assumed to have a discrete set of discontinuities. 

Particles reside at sites of Z, subject to the exclusion rule that admits at most one particle 
at each site. The dynamical rule is that a particle jumps from site i to site i + 1 at rate 
c(i/n) provided site i + 1 is vacant. Space and time are both scaled by the factor n and 
then we let n — > oo. We prove the almost sure vague convergence of the empirical measure 
to a density p(x,t), assuming that the initial particle configurations have a well-defined 
macroscopic density profile po. 

From known behavior of driven conservative particle systems a natural expectation would 
be that the macroscopic density p(x,t) of this discontinuous TASEP ought to be, in some 
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sense, the unique entropy solution of an initial value problem of the type 

(1.1) Pt + (c(x)f(p)) x = 0, p(x, 0) = p (x). 

Our proof of the hydrodynamic limit does not lead directly to this scalar conservation law. 
We can make the connection through some recent PDE theory in the special case of the 
two-phase flow where the speed function is piecewise constant with a single discontinuity. 
In this case the discontinuous TASEP chooses the unique entropy solution. We would 
naturally expect TASEP to choose the correct entropy solution in general, but we have not 
investigated the PDE side of things further to justify such a claim. 

The remainder of this introduction reviews briefly some relevant literature and then closes 
with an overview of the contents of the paper. The model and the results are presented in 
Section [2j 

Discontinuous scalar conservation laws. The study of scalar conservation laws 

(1.2) p t + F{x,p) x = 

whose flux F may admit discontinuities in x has taken off in the last decade. As with 
the multiple weak solutions of even the simplest spatially homogeneous case, a key issue is 
the identification of the unique physically relevant solution by means of a suitable entropy 
condition. (See Sect. 3.4 of [10] for textbook theory.) Several different entropy conditions 
for the discontinuous case have been proposed, motivated by particular physical problems. 
See for example [JJ El [U [9l [Til 02] • Adimurthi, Mishra and Gowda [3] discuss how 
different theories lead to different choices of relevant solution. An interesting phenomenon 
is that limits of vanishing higher order effects can lead to distinct choices (such as vanishing 
viscosity vs. vanishing capillarity). 

However, the model we study does not offer more than one choice. In our case the graphs 
of the different fluxes do not intersect as they are all multiples of f(p) = p(l — p). In such 
cases it is expected that all the entropy criteria single out the same solution (Remark 4.4 
on p. 811 of [2]). By appeal to the theory developed by Adimurthi and Gowda [T| we show 
that the discontinuous TASEP chooses entropy solutions of equation (II. ip in the case where 
c(x) takes two values separated by a single discontinuity 

Our approach to the hydrodynamic limit goes via the interface process whose limit is 
a Hamilton-Jacobi equation. Hamilton-Jacobi equations with discontinuous spatial depen- 
dence have been studied by Ostrov [12] via mollification. 

Hydrodynamic limits for spatially inhomogeneous, driven conservative parti- 
cle systems. Hydrodynamic limits for the case where the speed function possesses some 
degree of smoothness were proved over a decade ago by Covert and Rezakhanlou [8] and 
Bahadoran [5]. For the case where the speed function is continuous, a hydrodynamic limit 
was proven by Rezakhanlou in [13] by the method of [16] . 

The most relevant and interesting predecessor to our work is the study of Chen et al. 
[7]. They combine an existence proof of entropy solutions for (11. 2h under certain technical 
hypotheses on F with a hydrodynamic limit for an attractive zero-range process (ZRP) with 
discontinuous speed function. The hydrodynamic limit is proved through a compactness 
argument for approximate solutions that utilizes measure- valued solutions. The approach 
follows [SI [8] by establishing a microscopic entropy inequality which under the limit turns 
into a macroscopic entropy inequality. 
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The scope of [7] and our work are significantly different. Our flux F(x, p) = c(x)p(l — p) 
does not satisfy the hypotheses of [7|. Even with spatial inhomogeneities, a ZRP has 
product-form invariant distributions that can be readily written down and computed with. 
This is a key distinction in comparison with exclusion processes. The microscopic entropy 
inequality in [7] is derived by a coupling with a stationary process. 

Finally, let us emphasize the distinction between the present work and some other hy- 
drodynamic limits that feature spatial inhomogeneities. Random rates (as for example in 
[16| ) lead to homogenization (averaging) and the macroscopic flux does not depend on the 
spatial variable. Somewhat similar but still fundamentally different is TASEP with a slow 
bond. In this model jumps across bond (0, 1) occur at rate c < 1 while all other jump rates 
are 1. The deep question is whether the slow bond disturbs the hydrodynamic profile for 
all c < 1. V. Beffara, V. Sidoravicius and M. E. Vares have announced a resolution of this 
question in the affirmative. Then the hydrodynamic limit can be derived in the same way 
as in the main theorem of the present paper. The solution is not entirely explicit, however: 
one unknown constant remains that quantifies the effect of the slow bond (see [12]). [6] 
generalizes the hydrodynamic limit of [17] to a broad class of driven particle systems with 
a microscopic blockage. 

Organization of this paper. Section [2] contains the main results for the inhomogeneous 
corner growth model and TASEP. Sections [3] and 2] prove the limits. Section [5] outlines the 
explicit computation of density profiles for the two-phase TASEP. Section [6] discusses the 
connection with PDE theory. 

Notational conventions. N = {1,2,3, ... } and Z + = {0, 1,2,...}. The Exp(c) dis- 
tribution has density f(x) = ce~ cx for < x < oo. Two last passage time models appear 
in our paper: the corner growth model whose last-passage times are denoted by G, and 
the equivalent wedge growth model with last-passage times T. H(x) = 1j 0iO o)( x ) is the 
Heavyside function. C is a constant that may change from line to line. 



2. Results 



The corner growth model connected with TASEP has been a central object of study 
in this area since the seminal 1981 paper of Rost [14]. So let us begin with an explicit 
description of the limit shape for a two-phase corner growth model with a boundary along 
the diagonal. Put independent exponential random variables {l^j^eN 2 on the points of the 
lattice with distributions 



(2.1) 



Yt 



Exp(ci), 
Exp(ci A c 2 ) 
Exp(c 2 ), 



if i < j 
i£i = j 
if i > j. 



We assume that the rates satisfy c\ > c 2 . 
Define the last passage time 



(2.2) 



G(m, n) 



max 

7r£ll(m,n) 



Y v , (m,n)eN 2 
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Figure 1. Optimal macroscopic paths that give the last passage time con- 
stant described in Theorem 12.11 

where U(m, n) is the collection of weakly increasing nearest-neighbor paths in the rectangle 
[to] x [n] that start from (1,1) and go up to (to, n). That is, elements of U(m,n) are 
sequences {(1, 1) = v i, V2, ■ ■ • , v m+n -\ = (to, n)} such that Vi + \ — v j = (1, 0) or (0, 1). 



Theorem 2.1. Let the rates c\>C2> 0. Define c = c\/c2 > 1 andb = 2c—l — 2^c(c — 1) 
Then the a.s. limit 

&(x,y) = lim n" 1 G([nx\,[ny\) 



exists for all (x,y) G (0, oo) 2 and is given by 

ci 1 (V% + \fy) 2 > 



4c - (1 + b) 2 (1 + b) 2 - Acb 2 
x — tttt — h y- 



if < x < b 2 y 
if b 2 y < x < y 



Cl (l-b 2 ) y C!(l-6 2 ) 
S2 1 {V% + \/y) 2 > i/y<x<+oo. 
This theorem will be obtained as a side result of the development in Section [3l 



We turn to the general hydrodynamic limit. The variational description needs the fol- 
lowing ingredients. Define the wedge 

W = {{x, y) £ M 2 : y > 0, x > -y} 

and on W the last-passage function of homogeneous TASEP by 

(2.3) j(x, y ) = (V^+y + Vv) 2 - 

Let x(s) = (x\(s), X2(s)) denote a path in M 2 and set 

H(x,y) = {x G C([0, 1],W) : x is piecewise C\x(0) = (0,0), 

x(l) = (x,y), x (s) G W wherever the derivative is defined}. 
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The speed function c of our system is by assumption a positive lower semicontinuous 
function on M. We assume that at each i£l 

(2.4) c(x) = min < lim c(y), lim c(y) >. 

I y/-x y\x ) 

In particular we assume that the limits in (12, 4p exist. We also assume that c{x) has only 
finitely many discontinuities in any compact set, hence it is bounded away from in any 
compact set. 

For the hydrodynamic limit consider a sequence of exclusion processes rf 1 = (rjf(t) : i G 
Z, t G K+) indexed by n G N. These processes are constructed on a common probability 
space that supports the initial configurations {r] n (0)} and the Poisson clocks of each process. 
As always, the clocks of process rf 1 are independent of its initial state r] n (0). The joint 
distributions across the index n are immaterial, except for the assumed initial law of large 
numbers f|2. lOj) below. In the process rf 1 a particle at site i attempts a jump to i + 1 with 
rate c(i/n). Thus the generator of rf 1 is 

(2.5) L n f( V ) = c(xn" 1 )r ? (x)(l - v (x + l))(f(^ x+l ) - f( V )) 

for cylinder functions / on the state space {0, 1} Z . The usual notation is that particle 
configurations are denoted by r\ = (r](i) : i 6 Z) £ {0, 1} Z and 

{0 when i = x 
1 when i = x + 1 
rj(i) when i ^ x, x + 1 

is the configuration that results from moving a particle from x to x + 1. Let «/"(i) denote 
the number of particles that have made the jump from site % to site i + 1 in time interval 
[0,i] in the process r] n . 

An initial macroscopic profile po is a measurable function on R such that < po^) ^ 1 
for all real x, with antiderivative vo satisfying 

(2.6) v (0) = 0, v (&) - fo(a) = / po(x)dx. 



The macroscopic flux function of the constant rate 1 TASEP is 

p(l-p), if 0<p< 1 



(2.7) 

Its Legendre conjugate 




-oo, otherwise. 



inf {yp - f(p)} 

0<p<l 



represents the limit shape in the wedge. We orient our model so that growth in the wedge 
proceeds upward, and so we use g(y) = —f*(y)- It is explicitly given by 

-y, if y < -1 



(2.8) g(y) = sup {/(p) - yp} 

0<p<l 



1(1 -y)\ if - 1 < 2/ < 1 
0, ify>l. 
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For i£R define v(x, 0) = vq(x), and for t > 0, 



(2.9) v(x, t) = sup <^ v (w{0)) - / c(w{s)) g , ) ' ds 



w(-) I J 



* / w'(s) 



c(w(s)) 



where the supremum is taken over continuous piecewise C 1 paths w : [0, t] — > R that satisfy 
w(t) = x. The function v(x,t) is Lipschitz continuous jointly in (x,t). (see Section^]) and it 
has a derivative almost everywhere. The macroscopic density is defined by p(x, t) = v x (x, t). 

The initial distributions of the processes rf 1 are arbitrary subject to the condition that 
the following strong law of large numbers holds at time t = 0: for all real a < b 

1 l nb i ,6 

(2.10) lim - V t£(Q)= / po(x)dx a.s. 

n-too n ^— ' J a 
i=\na\+l 

The second theorem gives the hydrodynamic limit of current and particle density for 
TASEP with discontinuous jump rates. 

Theorem 2.2. Let c{x) be a lower semicontinuous positive function satisfying (|2.4p . with 
finitely many discontinuities in any compact set. Under assumption (|2.1U|) . these strong 
laws of large numbers hold at each t > 0: for all real numbers a < b 

(2.11) lim n _1 J7L i (nt) = vo(a) — v(a, t) a.s. 
and 

(2.12) lim- V ?7? (nt) = / p(x,t)dx a.s. 

i= [naj +1 

where v(x,t) is defined by (12.91) and p(x,t) = v x (x,t). 

Remark 2.3. In a totally asymmetric ET-exclusion with speed function c the state space 
would be {0, 1, ... , K} z with K particles allowed at each site, and one particle moved from 
site x to x + 1 at rate c{x/n) whenever such a move can be legitimately completed. Theorem 
12.21 can be proved for this process with the same method of proof. The definition of the 
limit (|2.9p would be the same, except that the explicit flux / and wedge shape g would be 
replaced by the unknown functions / and g whose existence was proved in [16j . 

To illustrate Theorem 12.21 we compute the macroscopic density profiles p(x,t) from con- 
stant initial conditions in the two-phase model with speed function 

(2.13) c(x) = c 1 (l-H(x)) + o 2 H(x) 

where H(x) = l[ 0jO o)( 2; ) is the Heavyside function and c\ > C2- (The case c\ < C2 can then 
be deduced from particle- hole duality.) The particles hit the region of lower speed as they 
pass the origin from left to right. Depending on the initial density p, we see the system 
adjust to this discontinuity in different ways to match the actual throughput of particles 
on either side of the origin. The maximal flux on the right is C2/4 which is realized on the 
left at densities p* and 1 — p* with 

p* = \ ~ \\fl-c 2 /a. 
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Corollary 2.4. Let c\ > c 2 and the speed function as in (|2.13|) . Then the macroscopic 
density profiles with initial conditions po(x,0) = p are given as follows. 

(i) Suppose < p < p*. Define r* = r*(p) = \- \ 4p(l - p)c 1 /c 2 - Then 

p if — oo < x < 

r* if < x < c 2 (l - 2r*)t 



(2.14) 



p(x,t) 



ih-^-J if C2 (l - 2r*)t < x < c 2 (l - 2p)t 
p if (1 — 2p)tc2 < x < +oo 




Figure 2. Density profile p(x,t) in the two-phase (ci > C2) TASEP when 
we start from constant initial configurations po(x) = p £ (0,p*). 



(ii) Suppose p* < p < \. Then 

P 

1 - p* 

(2.15) P(x,t) = l 1 



tc 2 



(hi) Suppose p > \. Define r* = r*{p) 



if — 00 < x < —tci(p — p*) 
if -t Cl (p- p*) < x < 

i/0 < x < (1- 2p)te 2 

i/ (1 — 2p)tc2 < x < +00 
i - — 4p(l — p)cijc\. Then 



(2.16) 



p if — 00 < x < —tc\{p — r*) 

1-r* i/ -ici(p-r*) < x < 
p if < x < +00 



Remark 2.5. Taking i — )■ 00 in the three cases of Corollary 12.41 gives a family of macroscopic 
profiles that are fixed by the time evolution. A natural question to investigate would be the 
existence and uniqueness of invariant distributions that correspond to these macroscopic 
profiles. 
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Figure 3. Density profile p(x,t) in the two-phase (ci > c 2 ) TASEP when 
we start from constant initial configurations po(x) = p G [p*, 5]. 





1-r* 




p 









T 

2 





c x {r* - p)t 

Figure 4. Density profile p(x,t) in the two-phase (ci > C2) TASEP when 
we start from constant initial configurations po{x) = p € (A, 1). 



Next we relate the density profiles picked by the discontinuous TASEP to entropy condi- 
tions for scalar conservation laws with discontinuous fluxes. The entropy conditions defined 
by Adimurthi and Gowda [lj are particularly suited to our needs. Their results give unique- 
ness of the solution for the scalar conservation law 



(2.17) 



j p t + (F(x,p)) x = 0, xeR,t>0 

\ p(x,o) = p (x), iei 

with distinct fluxes on the half- lines: 

(2.18) F(x,p) = H(x)f r (p) + (1 - H{x))h{p) 

where f T ,fe £ C 1 (1R) are strictly concave with superlinear decay to —00 as \x\ — > 00. 
A solution of (|2.17p means a weak solution, that is, p G L^ C (IR x M + ) such that for all 
continuously differentiable, compactly supported test functions £ C' (I x K+), 

r+OO r+OO / Q ± Q±\ r + OO 

(2.19) J J lp-^ + F{x,p)-£) dtdx + j p(x, 0)4>(x, 0) dx = 0. 
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(|2.19p is the weak formulation of the problem 

Pt + fr(p)x = 0, for x > 0, t > 

Pt + Mp)x = 0, forx<0,i>0 
f r (p(0+, t)) = fi(p(0-,t)) for a.e. t > 
p(ar,0) = po(x). 



(2.20) 



The entropy conditions used in [T] come in two sets and assume the existence of certain 
one-sided limits: 

(Ei) Interior entropy condition, or Lax-Oleinik entropy condition: 

(2.21) p(x+,t) > p(x-,t) for x / and for all t > 0. 

(Eb) Boundary entropy condition at x = 0: for almost every t, the limits p(0±,i) exist 
and one of the following holds: 



(2.22) 


fr(p(0+,t)) >0 and f e (p(0- 




(2.23) 


ti(p(0+,t)) <0 and /i(p(0- 


*)) < o, 


(2.24) 


/;(p(0+,t)) <0 and /;(p(0- 


t)) > o. 


Define 






G X ( P ) 


= i{ x > oyf*(p) + i{ x < o}/;(p) + i{x 


= 0} min 



where /* and /| are the convex duals of f r and fi. Set Vo(x) = f£ po(9) d9 and define 
(2.25) V(x,t) = S up{v (w(0))+ [ G w{s) {w'(s))ds) 

w(-) I JO ) 

where the supremum is over continuous, piecewise linear paths w : [0, t] — > M. with w(t) = 
x. 

Theorem 2.6. pQ Let p G L°°(M) and define V by (|2.25p . T/ien V is a uniformly Lipschitz 
continuous function and p(x, t) = V x (x, t) is the unique weak solution of (|2.2l)p that satisfies 
the entropy assumptions (Ei) and (Ef,) in the class L°°r)BVi oc and with discontinuities given 
by a discrete set of Lipschitz curves. 

It is easy to check that the two-phase density profile p(x, t) in Corollary 12.41 is a weak 
solution (in the sense of (]2.19p ) to the scalar conservation law (]2.17p with flux function 
F(x,p) = c(x)p(l — p). However we cannot immediately apply this theorem in our case 
since the two-phase flux function F(x,p) = (1 — H(x))cif(p) + H(x)c2f(p) is finite only 
for p £ [0, 1] and in particular is not C . We show how we can replace F(x, p) with F(x, p) 
in the above theorems in Section El In particular, we prove the following. 

Theorem 2.7. For p £ E define f r (p) = C2(l — p)p and fi(p) = c\(l — p)p to be the flux 
functions for the scalar conservation law (|2.20p . Let the initial macroscopic profile for the 
hydrodynamic limit be a measurable function < po(x) < 1. Then the macroscopic density 
profile p(x,t) from the hydrodynamic limit in Theorem \2.S\ is the unique solution described 
in Theorem W. 
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3. Wedge last passage time 

The strategy of the proof of the hydrodynamic limit is the one from [TB] and |17| . Instead 
of the particle process we work with the height process. The limit is first proved for the 
jam initial condition of TASEP (also called step initial condition) which for the height 
process is an initial wedge shape. This process can be equivalently represented by the 
wedge last-passage model. Subadditivity gives the limit. The general case then follows 
from an envelope property that also leads to the variational representation of the limiting 
height profile. In this section we treat the wedge case, and the next section puts it all 
together. 

Recall the notation and conventions introduced in the previous section. In particular, c(x) 
is a positive, lower semicontinuous speed function with only finitely many discontinuities 
in any compact set. Define a lattice analogue of the wedge W by 

(3.1) C = {(i,j)£Z 2 :j>l,i>-j + l} 

with boundary dC = {(i, 0) : i > 0} U {(i, -i) : i < 0}. 

For each n G N construct a last-passage growth model on C that represents the TASEP 
height function in the wedge. Let {jfj : G £}neN denote a sequence of independent 
collections of i.i.d. exponential rate 1 random variables. We need an extra index I to denote 
the shifting. Define weights 

-i 



(3.2) ^ = c (l_) , (i,j)e£ 



(' 



n 



For £ G Z and n G N assign to site (i,j) G C the random variable oj^jT^j. Given lat- 
tice points (a,b), (u,v) G C, Tl((a,b),(u,v)) is the set of lattice paths n = {(a,b) = 
(io,jo), (h,ji), ■••> (ipijp) = ( u , v )} whose admissible steps satisfy 

(3.3) (iiJi) - (ii-i,3i-i) G {(1,0), (0,1), (-1,1)}. 

In the case that (a,b) = (0,1) we simply denote this set by H(u,v). For (u, v) G £, £ G K 
and n G N denote the wedge last passage time 

(3.4) T^(u,v)= max £ 



mil \ ST^ n,t n 

' " ' max > bJ- '■ T; j 



with boundary conditions 

(3.5) T n ' e {u,v) = for (u,v) e dC. 

Admissible steps (|3.3p come from the properties of the TASEP height function. Notice 
that (0, 1) is in fact never used in a maximizing path. 

To describe macroscopic last passage times define, for (x, y) G W and g£l, 

(3.6) T«(x,y)= sup { t ; ds\. 

Theorem 3.1. For all q G M and (x,y) in the interior ofW 

(3.7) Mm n~ 1 T n ^([nx\ 1 [ny\)=T q {x,y) a.s. 
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Remark 3.2. In a constant rate c environment the wedge last passage limit is 
(3.8) lim -T n ([nx\,[ny\) = c-^y) = c- 1 + . 

n—too n 

The limit 7(2;, y) is concave, but this is not true in general for T°(x, y). In some special cases 
concavity still holds, such as if the function c(x) is nonincreasing if x < and nondecreasing 
if x > 0. 

To prove Theorem 13.11 we approximate c(x) with step functions. Let —00 = a\ < a 2 < 
... < cll-i < ciL = +00, and consider the lower semicontinuous step function 

L-1 L-1 

(3-9) c(x) = r m 1 (a m ,a m+1 )(x) + ^ nun{r m _i , r m }l {am} (x). 

m=l m=2 

Proposition 3.3. Let c(x) be given by (I3.9p . Then limit (|3.7p holds. 

On the way to Proposition 13.31 we state preliminary lemmas that will be used for pieces 
of paths. We write for the rate values instead of to be consistent with the notation in 
Theorem 12.11 

Lemma 3.4. Assume that there is a unique discontinuity 02 = for the speed function c(x) 
in (|3.9p . Then for y > 



lim n-'T^iO, [ny\) = — ^ -= -^-f ds a.s. 

mm{ci,c 2 } Jo c(0) 



n— >oo 



Proof. The upper bound in the limit is immediate from domination with constant rates 
<0). 

For the lower bound we spell out the details for the case c\ > 02- Let e > 0. To bound 
T n,0 (0, [ny\ ) from below force the path to go through points (0, 1), {( [rzi/ej ,(k — I) \nye\ ) : 
k = 1,... , le -1 ]} an d (0, [ny\). For 1 < k < [e _1 J let T n (R^) be the last passage time 
from ( \nye\ , (k — 1) \nye\ ) to ( [wyej , fc |_ n y £ J ) • R-k re fe rs to the parallelogram that contains 
all the admissible paths from ( \nye\ , (k — 1) \nye\ ) to ( [nyej , k [nye\ ). Each i?£ lies to the 
right of x = and therefore in the C2-rate area. (See Fig. [5J) 

Let < 5 < £C2 _1 7(0,y). A large deviation estimate (Theorem 4.1 in [15]) gives a 
constant C = C(c2,y,£,S) such that 

(3.10) P{2£(i#) < n^-SCO.y) - 5)} < e" c " 2 . 

By a Borel-Cantelli argument, for large n, 

r"' (0, \ny\) > T n (R n k )>n(le~ 1 \-l)(ec 2 ~ 1 1 (0,y)-5). 

k=l 

This suffices for the conclusion. □ 

Remark 3.5. This lemma shows why it is convenient to use a lower semi-continuous speed 
function. A path that starts and ends at the same discontinuity stays mostly in the low 
rate region to maximize its weight. This translates macroscopically to the formula for 
the limiting time constant obtained in the lemma, involving only the value of c at the 




discontinuity. If the speed function is not lower semi-continuous, we can state the same 
result using left and right limits. 

Lemma 3.6. Let a = < b < +00 be discontinuities for the step speed function c(x) and 
c{x) = r for a < x < b. Take z G [0,6]. Let T n ([nz\, [ny\) be the wedge last passage time 
from (0, 1) to ([nz\, \ny\) subject to the constraint that the path has to stay in the r-rate 
region (a,b) x (0,+oo), except possibly for the initial and final steps. Then 

(3.11) lim n- 1 T n (\nz\, [ny\) = r _1 7(z,y) a.s. 

n— >oo 

Same statement holds if b < z < 0. 

Proof. The upper bound lim n _1 T™( \nz\ , \ ny\) < r~ l ^(z,y) is immediate by putting con- 
stant rates r everywhere and dropping the restrictions on the path. For the lower bound 
adapt the steps of the proof in Lemma I3.4L □ 

Lemma [3.6l is a place where we cannot allow accumulation of discontinuities for the speed 
function. 

Before proceeding to the proof of Proposition 13.31 we make a simple but important ob- 
servation about the macroscopic paths x(s) = (x 1 (s), x 2 (s)), s G [0,1], in Ti(x,y) for the 
case where c(x) is a step function f|3. 9|) . 

Lemma 3.7. There exists a constant C = C(x,y,c(- — q)) such that the supremum in (13. 6p 
comes from paths inH(x, y) that consist of at most C line segments. Apart from the first and 



TASEP WITH DISCONTINUOUS JUMP RATES 



13 








Figure 6. A possible microscopic path described in Lemma 13.61 The path 
has to stay in the unshaded region. 

last segment, these segments can be of two types: segments that go from one discontinuity 
of c(- — q) to a neighboring discontinuity, and vertical segments along a discontinuity. 

Proof. Path x is a union of subpaths {xj} along which c(xj(s) — q) is constant, except 
possibly at the endpoints. Given such a subpath (xj(s) : tj < s < ij+i), concavity of 7 
and Jensen's inequality imply that the line segment <pj that connects x*(ii) to x.j(tj + i) 
dominates: 



Consequently we can restrict to paths that are unions of line segments. 

To bound the number of line segments, observe first that the number of segments that 
go from one discontinuity to a neighboring discontinuity is bounded. The reason is that 
the restriction x'(s) G W forces such a segment to increase at least one of the coordinates 
by the distance between the discontinuities. 

Additionally there can be subpaths that touch the same discontinuity more than once 
without touching a different discontinuity. Lower semi-continuity of c(-) and Jensen's in- 
equality show again that the vertical line segment that stays on the discontinuity dominates 
such a subpath. Consequently there can be at most one (vertical) line segment between 
two line segments that connect distinct discontinuities. □ 

Next a lemma about the continuity of T q . We write T q ((a, b), (x, y)) for the value in (|3,6p 
when the paths go from (a, b) to (x, y) G (a, b) + W. 

Lemma 3.8. Fix z,w > 0. Then there exists a constant C = C(z,w,c(- — q)) < 00 such 
that for all < 5 < 1 and < a < z 



Proof. Pick S G (0,1] and consider the point (z, 5) in W. For any x = (x 1 (s), x 2 (s)) G 
H(z, 6) set 





r9((a,0), (*,*)) -r*((a,0),(z,0)) < CV6, 



r«((-6, b), (-w, w + 6))- r«((-b, b), (-w, w 



)) < cVs. 



(3.14) 
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Let e > and assume that 4> = 4> 2 ) £ "^(-z, is a path such that T q (z, 5) — I(<p, q) < e. 
Lemma 13.71 implies that we can decompose <f> into disjoint linear segments 4>j so that <f> = 
Ylj=i 4>j an d <Pj '■ [ s j-i> s j] W. Here Ylj 4>j means path concatenation. 
We can find segments <j>j{h\, 1 < & < iV, such that 

^}(fc)( s i(fc)-i) < ^i(fc)( s i(fc))' 0j{k)( s j{k)) = 0i(ik+i)( s i(fc+i)-i)> 

= 0, and 4>^N^( s j(N)) = z - ^ n °t ner words, the projections of the segments 
4>j(k) cover the interval [0, z] without overlap and without backtracking. 

We bound the contribution of the remaining path segments to I(4>,q). Let J = 
\Jk=i i s j{k)i s i(fc+i)-i] be the leftover portion of the time interval [0,1]. The subpath 
4>(s), s E [sj(k)j s j(k+i)-x]i (possibly) eliminated from between 4>j(k) an d <f>jfk+l) satisfies 
^(sj^) = ^(^(fc+i)-!)- Note tnat 7( a ' b) <2a + 46 for (a, i))eW and J 1 ( ( /) 2 ) / (s) a?s = 5. 
We can bound as follows: 

j c(4>Hs)-q) dS - C Jj ^ ] (S) ' { ^ } ^ 

<C^(2(0 1 )'( S )+4(0 2 )'( S ))& 

<C f 2(0 1 )'(s)ds + C / 4(0 2 )'(s) ( is 

Jo 

(3.15) = + 4C5. 

Set tfc = Sj(fc)_i < nfc = Sj(fc). Define a horizontal path w from (0,0) to (2,0) with 
segments 

(3.16) w k (s) = (0] (fe) (s),O), for t k <s<u k , 

and constant on the complementary time set J. 
To get the lemma, we estimate 

ij c(0 A (s) - g) 7 [0,1] V c (0 ( s ) - 

N 



<cs + J2 (^-(k), ?) - 9)) + r"(z, o) 

fe=l 

< C( 5 + c"^ / ( 7 (^ (fc) ( a ))- 7 (u4( fl )))cfa + r9(z,o) 

J„ 1 -'tfc 



fc=l Jt >= 

<C8 + C>Y,( / (0% k )( s )ds+ 



k=l 



< C5 + C'j2<yf^^ 2 )' m (s)ds^ 2 ( jT + (0 2 ); (fc) (.)) da) 2 0) 
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(0,0) (*,0) (z + 5,0) 

Figure 7. A possible macroscopic path from (0,0) to (z,5). The dotted 
vertical lines are discontinuity columns of c(— q). The error from eliminating 
the segments outside the two vertical dashed lines and from eliminating 
pathologies (like the circled part) is of order 5 and a comparison with the 
horizontal path leads to an error of order \f5. 



E / W*) x 

fc=i Jt *> ' 

E / (^y m (s) + (^ 2 y m (s)) ds) +H(z, ) 

1 1 Jtk / 



fc=l ^ tfe 

<^+c^^ / ((<A 1 );. (fc) ( S ) + (<A 2 ); (fc) ( S ))^+H(z,o) 
fe=i 

< C5 + c'^V^T^ + r q {z, o) 

< C<5 + C7'^Vi + r 9 (z,0). 

The first inequality (|3.12p follows for a = by letting e go to 0. It also follows for all 
o G [0, z] by shifting the origin to a which replaces z with z — a. 

For the second inequality f|3. 13|) the arguments are analogous, so we omit them. □ 



Corollary 3.9. Fix (x,y) £ W. Then there exists C = C(x,y,c(- — q)) < oo such that 
for all0<5<l 

(3.17) T q (x,y + 5)-T q (x,y) <CV6. 



Proof. Let A((a, b), (x, y)) be the parallelogram with sides parallel to the boundaries of the 
wedge, north-east corner the point (x,y) and south-west corner at (a, b). If (a, b) = (0,0) 
we simply write A(x,y). 

Let e > 0. Let 4> £ a path such that T q (x, y + 5) — I{4> £ ,q) < e. Let u be the point where 
(j> £ first intersects the north or the east boundary of A(x,y). Without loss of generality 
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assume it is the north boundary and soti= (a, y) for some a G [— y, x]. Then, 

r«(x,y + 6)-e<I(<f> £ ,q) 

<T«(a,y) + ri((a,y),(x,y + 5)) 

= r«(o, y) + H((a, y), (x, y)) + H((a, y), (x, y + 6)) - H((a, y), (x, y)) 

(3.18) < r«(x, y) + H((a, y), (x, y + 5)) - T«((a, y), (x, y)). 

The last inequality gives 

(3.19) r*(z, y + 5) - H(x, y) < H((a, y), (x, y + 5)) - H((a, y), (x, y))+s<CVd + s 
by Lemma 13.81 Let e decrease to to prove the Corollary. □ 

Proof of Proposition HOI Fix (x,y) in the interior of W. For x = (ar(s), x 2 (s)) G H(x,y) 
set 



(3.20) 



We prove first 



(3.21) Km ^^'^(LnxJ, [nyj) >T«(x,y) = sup J(x,g). 

«->oo x(-)e'H(x,y) 

It suffices to consider macroscopic paths of the type 

H 



(3.22) x( S ) = ^ Xj ( S )l[ 



Sj,S j+1 ) 



S 



3=1 



where H G N, Xj is the straight line segment from x(sj) to x(s J+ i), c(xi(s) — g) = r m . is 
constant for s G (sj,s J+ i), and by continuity Xj(sj + i) = Xj +1 (s J+ i). 

Let 7r n be the microscopic path through points (0,1), { |_ nx j(sj)J : 1 < J < K} and 
( [nxj , [ny\ ) constructed so that its segments ir™ satisfy these requirements: 

(i) 7T™ lies inside the region where = is constant, except possibly for the initial 
and final step; 

(ii) 7r™ maximizes passage time between its endpoints [nx.j(sj)\ and \nxj+i(sj+i)\ sub- 
ject to the above requirement. 

Let 

(3.23) T^=msx £ ^r i;k 

denote the last-passage time of a segment subject to these constraints. Observe that the 
proofs of Lemmas 13.41 and 13.61 do not depend on the shift parameter q, therefore 

hm n-WM = 7(x^)-x J+1 ( Sj)) = p 7^ («)) ^ 
3 r mj J Sj c(xj(s)-g) 
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Adding up the segments gives the lower bound: 




Km n~ l T n '\- n ^{[nx\, [ny\) > lim ^ n -^T^ m] 




p+i 7(x;-(s)) 
J s . c(x 1 (s)-q) 



ds = /(x, q). 



Now for the complementary upper bound 



(3.24) lim n- l T n '\- nq ^{[nx\ : [ny\) < T q (x,y). 



Each microscopic path to ( \nx\ , [ny\ ) is contained in nA for a fixed parallelogram A C W 
with sides parallel to the wedge boundaries. Pick e > 0. Let r, > be a lower bound on 
all the rate values that appear in the set A. Find 5 > such that ^(v) — 7(1^)! < er* for 
all v,w £ A with \v — w\ < 5 and 5 < 1 so that Corollary 13.91 is valid. 

Consider an arbitrary microscopic path from (0,1) to ([^J, |_n2/J ) . Given the speed 
function and q, there is a fixed upper bound Q = Q(x,y) on the number of segments of the 
path that start at one discontinuity column ([raajj + [nq\) x N and end at a neighboring 
discontinuity column ([naj-tij + [nq\) x N. The reason is that there is an order n lower 
bound on the number of lattice steps it takes to travel between distinct discontinuities in 
nA. 

Fix K 6 N and partition the interval [0, y] evenly by bj = jy/K, < j < K, so that 
y/K < 5/Q. Make the partition finer by adding the y— coordinates of the intersection points 
of discontinuity lines {en + q} x M + with the boundary of A. 

Let 7r n be the maximizing microscopic path. We decompose ir n into path segments 
{7r™ : < j < M n } by looking at visits to discontinuity columns (L^ctjJ + \nq\) x N, both 
repeated visits to the same column and visits to a column different from the previous one. 
Let {0 = bk < bk x < bk 2 < ... < t>k Mn _ 1 < ^k Mn = y} be a sequence of partition points and 
{0 = xq, x\ = a mi +q, X2 = a m , 2 +q, . . . , xm„ = x} & sequence where Xj for < j < M n are 
discontinuity points of the shifted speed function c(- — q). We can create the path segments 
and these sequences with the property that segment 7r™ starts at ([nxj\,l) with / in the 
range [nb k .\ < I < [nb k . +1 \ and ends at ([nxj + i\, I') with [nb kj+1 \ < V < [nb kj+1+1 \. In an 
extreme case the entire path -zr™ can be a single segment that does not touch discontinuity 
columns. 

In order to have a fixed upper bound on the total number M n of segments, uniformly in 
n, we insist that for < j < M n — 1 the labels satisfy: 

(i) For odd j, tt™ starts and ends at the same discontinuity column (L^XjJ , • ). The rate 
relevant for segment 7r" is = c{a mj ) . 

(ii) For even j, ir™ starts and ends at different neighboring discontinuity columns, and 
except for the initial and final points, does not touch any discontinuity column and visits 
only points that are in a region of constant rate m ■ . 

The above conditions may create empty segments. This is not harmful. Replace Q with 
2Q + 2 to continue having the uniform upper bound M n < Q. 

Let T{tt^) be the total weight of segment 7r™. Let 7r" be the maximal path from 

(|/JXjJ, \ nb kj \) to (|/JXj + iJ, [nb k]+1+ i\ ) in an environment with constant weights Uij = rj l 
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everywhere on the lattice, with total weight Tj 1 . Tj 1 > T(7r"), up to an error from the end- 
points Of 7T™. 



Theorem 4.2 in [15] gives a large deviation bound for T™. Consider a constant rate r 
environment and the maximal weight T((\nui\, [n^ij), ([mt2.|> L 7 ^])) between two points 
(u\,vi) and (^2,^2) such that their lattice versions can be connected by admissible paths 
for all n. Then there exists a positive constant C such that for n large enough, 



< e~ Cn . 



(3.25) P^T(([nui\,[nvi\),([nu2\,[nv2\)) > nr lr y(u2 — ui, v% — vi) + ne| 

There is a fixed finite collection out of which we pick the pairs {(xj, 6^.), (zcj+i, &fc 3+1 +i)} 
that determine the segments fr™. By (|3.25p and the Borel-Cantelli lemma, a.s. for large 
enough n, 

(3.26) T/ < nr- 1 7 (^+i - sj, ^ J+1 +i - b k .) + ne for < j < M n . 



Define 5 X > by j, + 5 X = E^V^^+i+i " Since y = E^V^^+i " and 

by the choice of the mesh of the partition we have 5\ < M n 5/Q < 5. Think of 

(xj + i — Xj, 6fc. +1 -fi — 6t ), < j < M n , as the successive segments of a macroscopic path 
from (0,0) to (x, y + <5i). 

For sufficiently large n so that (|3.26p is in effect, 

T n ' ^ ( [nx\ , [ny\ ) < ^ T™ < n ^ ^(sj+l " ^ , &ft J+1 +l " &*, ) + 

< nT q (x,y + 81) + nQe 

< nT q (x, y) + nC\Tb~ + nQe. 

The last inequality came from Corollary 13.91 Let 6 — > 0. Since e was arbitrary the upper 
bound ([351]) holds. □ 

Proof of Theorem \3.1[ Fix (x,y). For each e > we can find lower semicontinuous step 
functions c\ and C2 such that ||ci — C2II00 < e and on some compact interval, large enough 
to contain all the rates that can potentially influence T q (x,y), c\{x) < c(x) < 02(2). When 
the weights in (|3.2|) come from speed function Cj let us write Tj for last passage times and 
Ti for their limits. An obvious coupling using common exponential variables {Tij} gives 

T^ nq \\nx\, \ny\) > T n ^ (\nx\ , [ny\ ) > ([nx\ , [ny\ ). 

Letting a > denote a lower bound for c(x) in the compact interval relevant for (x,y), we 
have this bound for x € %{x,y): 

J lci(xi(s) - c 2 (xi(s)-g)J 7 cf(xi(s)-g) 

< ea~ 2 j(x,y). 

Therefore the limits also have the bound 

0<Tl(x,y)-Tl(x,y) <C(x,y)e. 

From these approximations and the limits for Tj in Proposition 13.31 we can deduce Theorem 
l3Tl □ 
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Proof of Theorem \2.1\ We can construct the last passage times G(x, y) of the corner growth 
model (12. 2p with the same ingredients as the wedge last passage times T n '°(x, y) of fl33J, by 
taking Y^j) = j T i-j j- Then T n,0 (x, y) = G(x + y, y) and we can transfer the problem 
to the wedge. The correct speed function to use is now c(x) = c\\{x < 0} + C2l{x > 0}. 
In this case the limit in Theorem 13.11 can be solved explicitly with calculus. We omit the 
details. □ 

4. Hydrodynamic limit 

In this section we sketch the proof of the main result Theorem 12.21 This argument is 
from [III Q2]. 

4.1. Construction of the process and the variational coupling. For each n G N we 
construct a Z- valued height process z n (t) = (zf(t) : i £ Z). The height values obey the 
constraint 

(4.1) 0<zf +1 (t)-zr(t)<l. 

Let {£>"} be a collection of mutually independent (in i and n) Poisson processes with rates 
cf given by 

(4.2) c? = c(n-U), 

where c(x) is the lower semicontinuous speed function. Dynamically, for each n and i, the 
height value z™ is decreased by 1 at event times of X>", provided the new configuration does 
not violate (14. ip . 

After we construct z n (t), we can define the exclusion process r] n (t) by 
(4-3) = *?(*) -,?_!(*). 

A decrease in zf is associated with an exclusion particle jump from site i to i + 1. Thus 
the z n process keeps track of the current of the ry^-process, precisely speaking 

(4.4) J?{t) = <(0) - *?(*)■ 

Assume that the processes z n have been constructed on a probability space that supports 
the initial configurations z n (0) = {zf (0)) and the Poisson processes {T>f} that are indepen- 
dent of (zf(0)). Next we state the envelope property that is the key tool for the proof of 
the hydrodynamic limit. Define a family of auxiliary height processes {£ n ' fc : n G N, k G Z} 
that grow upward from wedge-shaped initial conditions 



(4.5) C' fc (0) 



0, if i > 
-i, if i < 0. 



The dynamical rule for the £ n ' k process is that jumps up by 1 at the event times of 
T>2 +k provided the inequalities 

( 4 -6) C t ' <fili and <£4i + l 

are not violated. In particular ^ n ' fc attempts a jump at rate c™ +fe . 
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Lemma 4.1 (Envelope Property). For each n G N, for all i G Z and t > 0, 
(4.7) z?(t) = supK(0)-ei.(i)} a.*. 

Equation (|4.7p holds by construction at time t = 0, and it is proved by induction on 
jumps. For details see Lemma 4.2 in [16] . 



4.2. The limit for £. For q, x G R, t > and for the speed function c(x), define 

(4.8) g q {x, t) = inf {y : (x, y) G W, r 9 (x, y) > t} . 

T q (x, y) defined by (j3.6H represents the macroscopic time it takes a £-type interface process 
to reach point (x, y). The level curve of T q given by g g (-,t) represents the limiting interface 
of a certain ^-process, as stated in the next proposition. 

Proposition 4.2. For all q,x G M. and t > 

(4.9) lim n~ l & (nt) = g~ q (x, t) a.s. 

Recall the lattice wedge L defined by (|3.ip . For (i, j) G £ U dC, let 

(4.10) L n ' k (i,j)=mi{t>0:^' h (t)>j} 

denote the time when £ n ' fc reaches level j. The rules (|4.5p ~ (|4.6p give the boundary conditions 



(4.11) L n > k (i,j)=0 for (i,j)€d£ 
and for G C the recurrence 

(4.12) L n ' k (i,j) = max{L"' fc (i - l,j),L n > k (i,j - l),L"- fc (i + l,j - 1)} + ^ 

where j3™'- is an exponential waiting time, independent of everything else. It represents the 
time £™ ,k waits to jump, after and its neighbors Ci~-v have reached positions that 
permit £™' k to jump from j — 1 to j. The dynamical rule that governs the jumps of 
implies that f3™' k has rate c™ +fc . 

Equations (13. 4h . (13. 5h . (14. lip , and (14.121) . together with the strong Markov property, 
imply that 

(4.13) {L n > k (i,j) : (i,j) G £ U <9£} = {T n >~ k (i,j) : eCUdC}. 

Consequently Theorem 13.11 gives the a.s. convergence n - l L n '\- n ^{\nx\, [ny\) -> T-*(x,y), 
and this passage time limit gives limit (|4.9p . 

Proof of Theorem \2.2l Given the initial configurations rj n (0) = {rjf(0) : i G Z} that appear 
in hypothesis (12.101) . define initial configurations z n (0) = {z™(0) : % G Z} so that 2q (0) = 
so that (|4.3p holds at time i = 0. Hypothesis (I2.10p implies that 

(4.14) lim n~ x zl ■ , = v {q) a.s. 
for all gGM, with t;q defined by (12.61) . 



TASEP WITH DISCONTINUOUS JUMP RATES 21 

Construct the height processes z 11 and define the exclusion processes rf 1 by (I4.3p . Define 
v(x,t) by (|2.9p . From (|4.3p ~ (14.4[) we see that Theorem 12.21 follows from proving that for all 

x £R,t£R + , 

(4.15) \\m n~ l z?^Ant) = v(x,t) a.s. 

n-foo L"^- 1 

Rewrite (14. 7p with the correct scaling: 



(4-16) n- 1 ^ (nt) = sup {n" 1 ^ (0) - n" 1 ^.^ (n*)} . 

The proof of (|4.15jl is now to show that the right-hand side of (|4.16p converges to the 
right-hand side of (j2.9|) . 

From (|4.14p . ()4.16p and (14. 9h we can prove that a.s. 

(4.17) lim n- 1 z? nxl (nt) = snp{v (q)-g- q (x-q,t)}=v(x,t). 



The argument is the same as the one from equations (6.4)-(6.15) in [16] so we will not 
repeat it here. 

Using (|3.6p and f|4.8j) we can rewrite v(x,t) as 

(4.18) v(x,t) = sup < t>o(<?) — y '■ 3x € %(x — q, y) such that / — — — -ds>t>. 

q , y m L " Jo c(xi(a) + g) J 

The final step is to prove f (x, t) = w(x, t). The argument is identical to the one used to prove 
Proposition 4.3 in [T7] so we omit it. With this we can consider Theorem 12.21 proved. □ 

5. Density profiles in two-phase TASEP 

This section proves Corollary 12.41 assuming c(x) = (1 — H(x))c\ + H(x)c2, c\ > C2 and 
Po(x) = p € (0, 1), we use variational formula (|2,9p to obtain explicit hydrodynamic limits. 

Remark 5.1. In light of Corollary [221 one can (instead of doing the following computations) 
guess the candidate solution for the scalar conservation law (|2.20p and then check that it 
verifies the entropy conditions (|2.22p - (|2.24p . The following computations do not require 
any knowledge of p.d.e. theory or familiarity with interface problems so we present them 
independently in this section. 

Let 

(5.1) C (x, t, q) = {w 6 C([0, t],M) : w piecewise linear, w(0) = q, w(t) = x} . 

To optimize in (|2.9p we use a couple different approaches for different cases. We outline 
this and omit the details. 

One approach is to separate the choice of the starting point q of the path. By setting 

(5.2) I(x,t,q)= inf ( I c(w(s))g ( W [ S \^ \ ds\ 



(12. 9p becomes 

(5.3) v(x, t) = sup {v (q) - I(x, t, q)} 
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We distinguish four cases according to the signs of x, q. Set 

(5.4) R+{x, t) = sup {vo(q) — I(x, t, q)} , if x > 0, 

(5.5) L-(x, t) = sup {vo(q) — /(x, t, q)} , if x < 0. 

These functions are going to be used in Cases 1 and 2 below {qx > 0) where we can compute 
I(x,t,q) directly. 

However, there are values (x, t, q) for which the (/-derivative of the expression in braces in 
(15. 3p is a rational function with a quartic polynomial in the numerator. While an explicit 
formula for roots of a quartic exists, the solution is not attractive and it is not clear how 
to pick the right root. Instead we turn the problem into a two-dimensional maximization 
problem. 

If qx < the optimizing path w crosses the origin: w(u) = for some u. It turns out 
convenient to find the optimal q for each crossing time u. For Case 3 (q < 0, x > 0) set 

(5.6) ®(u,q) = qp- ciug ( —) - c 2 {t -u)g 



UC\ J \ (t — U)C2 

and 

(5.7) L + (x,t) = sup <&(u,q). 

q<0,ue[0,t] 

For Case 4 (q > 0, x < 0) the obvious modifications are 

/ — q \ fx 

(5.8) ^(u,q) = qp- c 2 ug [ -ci(i-ti)-' 



UC2 J \ (t — U)C\ 



and 

R-{x,t) = sup ^(u,q). 

g>0,ue[0,t] 



Rewrite (|5.3p using functions R±, L±: 

(5.9) v(x,t) = max{i?+(x,t), L + (x,t)}l{x > 0} + max{i?_ (x, i), L_ (x, t)}l{x < 0}. 

Proof of Corollary \2.4\ We compute the functions R±, L±. The density profiles p(x,t) are 
given then by the x-derivative of v(x,t). 

Case 1: x > 0, q > 0. Since c% < ci, the minimizing w of I(x,t,q) is the straight line 
connecting (0,(7) to (t,x). In particular, 

(5.10) I(x,t,q) = c 2 tg 



tC2 

Then the resulting R + (x,t) is given by 

'-tc 2 g(j^) ifp<3, x<tc 2 (l-2 / o) 

/9X — tC 2 p(l — P), if p < |, X > tC2(l — 2p) 

^px-tc 2 p{\- p), ifp>|, 



(5.11) i?+(x,t)=< 
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Case 2: x < 0, q < 0. The minimizing path w can either be a straight line from (0, q) 
to (t, x) or a piecewise linear path such that the set {t : w(t) = 0} has positive Lebesgue 
measure. This last statement just says that the path might want to take advantage of the 
low rate at x = 0. We leave the calculus details to the reader and record the resulting 
minimum value of I(x, t, q). Set B = yc\{ci — c%). 

' -9 c i (-[ _ B\ 2 , ( f _ \x\-g \ c 2 _ xci (i , B \ 2 
4B \ L ci y y B J 4 4B \ L ^ Cl ) ' 



(5.12) J(a;,t, g ) = 4 



( ^ 



when - - -v/N") 2 < 9; < x < 

otherwise 



The corresponding function L_(x,i) is given by 
(5.13) 

px - tc\p{\ - p), 



L-(x,t) = < 



px - tc\p(\ - p), 

-(* + §)? + ft( 1 + f) 

pX - tCip(l - p), 

-(* + §) f + 
-(* + §)? + + 

px - citp(l - p), 



ci y ' 2 

2 



< p < p*,x e r 

P*<P<^, x < -tci{p - p*) 
P*<P<\, x>-tc 1 {p-p*) 

t; < P < 1 — p* , X < —tC\(p — p*) 

1 <p<l-p*, -tc 1 (p-p*)<x 

1 - p* < p < 1, -St < X 

l-p*<p<l, -cit(2p - 1) < x < -Bt 
1 - p* < p < 1, x< -cit(2p - 1) 



Case 3: x > 0, q < 0. Abbreviate D = c| — 4ciC2p(l — p). First compute the ^-derivative 

l 2_ 

2 2mci 



(5.14) 



q < —uc\. 



If p > 1/2 then & q is positive and the maximum value is when q = so we are reduced to 
Case 1. If p < 1/2 the maximizing g = uci (2p — 1) . Then 

F(u) = 2uci(p - ±)) = -ncip(l - p) - c 2 (t - u)g { jj-^T^ j , 



with u-derivative 



dF . . c 2 / 

- = - Cl p(l - p) + T ^1 - 2; . 

Again we need to split two cases. If p < p* (equivalently D > 0) and x < tv^D, the 
maximizing u = t — x/y/D, otherwise u = 0. If p* < p < ^ the derivative is negative so the 
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maximizing u is still u = 0. Together, 



' -tc lP {\ - p) + x{\ - ^ ), P <p*,x<tVD 



(5.15) 



L+(:r,i) 



p < p*,x > t\f~D 
p*<p<l. 



Case 4: x < 0, q > 0. We treat this case in exactly the same way as Case 3, so we omit 
the details. Here we need the quantity D\ 

(5.16) R-(x,t) = 



(ci) — 4ciC2p(l — p) and we compute 
P<\ 



tc 2 p(l- p)+x(± + ^), \< p,-ty/Dl<X 



.- tc ^(j^), \<P<l,x<-t^fD[ 
Now compute v(x,t) from (|5.9p . We leave the remaining details to the reader. 



□ 



6. Entropy solutions of the discontinuous conservation law 

For this section, c(x) = (1 — H(x))c\ + H(x)c2, h(p) = p(l — p) and set F(x,p) = c(x)h(p) 
for the flux function of the scalar conservation law (|2,17p and F(x,p) = c(x)f(p) for the 
flux function of the particle system, where / is given by (|2.7p . (The difference between F 
and F is that the latter is — oo outside < p < 1.) 

In p] the authors prove that there exists a solution to the corresponding Hamilton-Jacobi 
equation 

' V t + cxh(V x ) = 0, if x < 0,t > 

(6.1) I V t + c 2 h(V x ) = 0, ifz>0,£>0 

, V(x,0) = Vo(x) 

such that solves the scalar conservation law (I2.17P with flux function F(x, p) and V x 
satisfies the entropy assumptions (Ei), (E^). V(x,t) is given by 

(6.2) V(x,t) =supJV (™(0)) + / (c(w(s))h)*{w'{s))ds 

w(-) { Jo 

where the supremum is taken over piecewise linear paths w G C([0, t], M) that satisfy w(t) = 
x. 

To apply the results of pQ to the profile p(x,t) coming from our hydrodynamic limit, 
we only need to show that the variational descriptions match, in other words that we can 
replace F with F and the solution is still the same. 



Proof of Theorem \2.7\ Convex duality gives (c(x)/)* (y) = c(x)f* (y/c(x)) and so we can 
rewrite (|2,9h as 



•3) 



v(x,t) = sup < ^o(^(0)) 



(c(w(s))fy(w'(s))ds 
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Observe that for all y 6 R 

(6.4) (c(x)fr(y)>(c(x)h)*(y), 

with equality if and only if y € [— ci,C2] Since the supremum in ()6.2f) and (16.3P is taken 
over the same set of paths, (|6.4p implies that 

(6.5) V(x,t) <v(x,t). 

The proof of the theorem is now reduced to proving that the supremum in (|6.3p is achieved 
when w'(s)c(w(s))~ 1 £ [—1,1], giving V(x, t) =v(x,t). 
To this end we rewrite v(x,t) once more, this time as 

v(x,t) = max{i? + (x, t), L + (x, t)}l{x > 0} + max{i?_(x,i),L_(x,t)}l{x < 0} 

where the functions R± , L± (as in the proof of Corollary 12. 4p are defined by 

(6.6) R+(x, t) = sup {vo(q) — I(x, t, q)} , if x > 0, 

(6.7) L-(x,t) = sup{vo(q) — I(x,t,q)} , if x < 0, 

g<0 

where I(x,t,q) is as in ()5.2|) . and 



(6.8) L+(x,t) = sup <v (q)-ciug( ) - c 2 (t - u)g r— ^ if x > 0, 

9 <0,u6[0,t] I Vuci/ \{t-U)C 2 ' 

and 

(6.9) R_(x,t) = sup i « (g) ~ C2Uff ( ) -ci(t- u)g I T ) }>, x < 0. 

«>o,ue[o,t] I V uc 2/ V(* — M ) c i . 

It suffices to show that the suprema that define R± , L± are achieved when 

(6.10) u/(s)c(w(s)) _1 € [-1,1]. 

We show this for L+. The remaining cases are similar. In (|6.8p . as before, it is the time for 
which w(u) = 0. Let $(u, q) denote the expression in braces in (|6.8j) with g-derivative 

Po{q) - I - 2^T> _,uc i < <? < 

Observe that if $> q (u,q) = for some q* = q*(u) then also q* maximizes $. Otherwise 
the maximum is achieved at and we are reduced to a different case. Assume that q* exists. 
Then by (IBTTTI) 

* 

(6.12) ^- = (l-2p (g*))ci <ci. 

n 

Therefore, the slope of the first segment of the maximizing path w satisfies (|6.1(JI) . 

The slope of the second segment is x[t — u) _1 . Assume that the piecewise linear path 
w defined by u and q* is the one that achieves the supremum. Also assume u > t — xc^ 1 - 
Consider the path w with w(0) = q* and w(t — xc^ 1 ) = 0. Since g is decreasing, we only 
increase the value of Hence the supremum that gives L + cannot be achieved on w and 
this gives the desired contradiction. □ 



(6.11) %(u, 
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